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Abstract 

In the low energy region chiral perturbation theory including virtual photons 
^' is used to derive the structure of the generating functional. The work we 

^ . do is performed within the three flavor framework and reaches up to next- 

to-leading order. An Euclidian cut-off is introduced to separate the long- 

'k>( \ and short-distance contributions. The long-distance part is evaluated in the 

\-{ ' ChPT framework up to 0{p^). The short-distance part is achieved partly 

from the work of J. Bijnens and J. Prades P, through perturbative QCD 
and factorization. A matching is made and the finite parts of the low-energy 
constants (LECs) is determined by using already existing data of the QCD 
LECs up to order / P, 0, 0- 
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1 Introduction 

Masses of hadrons have always been a problem to calculate throughout the history of 
modern particle physics. The problem lies in the fact that the particles are believed to 
have three different mass contributions; the free quark mass from the Higgs mechanism, 
the constituent mass and perturbative interactions with 7, W and Z. The last contribution 
is the one that we are concentrating on. These are called radiative corrections and involve 
interactions from gauge bosons. For the light mesons that we are going to investigate 
it is enough to use photon corrections, in an effective low energy model. This model is 
derived from the Standard Model (SM) which is a relativistic quantum field theory that 
give rise to the strong interaction described by Quantum Chromo Dynamics (QCD) and 
the electroweak interaction described by Quantum Flavor Dynamics (QFD). QCD works 
very well and agrees with experiments at high energies and short distances. It is much more 
difficult to examine QCD at low energies and long distances. One of the most successful 
approaches is based on the fact that QCD has a SU{nF)L x SUinp)^ chiral symmetry for 
Up massless flavors. To get an effective field theory model for low energy, a perturbation 
of the symmetry above is performed. The perturbation is called chiral perturbation theory 
(ChPT). This theory can be described by certain fields which correspond to the low energy 
spectrum of light pseudoscalar mesons, the tt, K, and rj octet. When Cem is added to Cqcd 
the effective Lagrangian for QCD and Electromagnetic (EM) interactions using Weinberg's 
power counting scheme jl| can be written as: 

10 2 14 

CeffecUve = £2 + A + " " " = ^2(0, Fq, Bq) + J^ L^Ol + J] H.Ol + J^ K^^l + ■■■ (1) 

j=l i=l j=l 

where the subscript has been ordered in momenta, mass and electric charge. The O's 
represent the term in the effective Lagrangian that are connected to the i:th LEC (see 
Eq. fl44|) ) . The lowest order Lagrangian is given by the non-linear sigma model coupled to 
external fields [S]|Hj. It contains three free parameters in the chiral limit: the scalar quark 
condensate Bq, the pion-decay constant Fq and the lowest order EM coupling constant C. 
The leading order expansion 0{p'^) implies tree level diagrams with C2 vertices. The next- 
to-leading order expansion is used to derive the electromagnetic LECs Ki by computing 
photon loops from the leading order Lagrangian £2 and adding contributions from the 
next-to-leading order Lagrangian £4. Doing the loop calculations will involve all the 29 
LECs^ in Eq. (P) that are connected to the couplings of different vertices in the diagrams, 
more about this in section 2.5. This is not a problem to us since the purely QCD part 
connected to Bq, Fq, Li and Hi, have been determined experimentally or theoretically. The 
constants of order 0{p'^) contain a divergent part and a constant finite part. The divergent 
part can be isolated with the help of dimensional regularization and MS regularization. 
For the K^ part this was recently done by Urech in 0, but the finite part of the Ki is 
almost entirely unknown. The goal of this thesis is to determine the K^ as a function of 
the Li and Hi. 



^ There exist a few more constants up to order 0{p'^), but these are of no importance to us 



2 Theory 

One of the greatest achievements in physics during the 20th century is the Standard model 
of particle physics. It revolutionized the view of elementary particles. Protons and neutrons 
that build up atoms (and electrons, of course) are not considered to be the fundamental 
particles any more. Instead they consist of point-like particles called quarks. Beside the 
quarks exist fundamental particles called leptons. In total there are six quarks (u, d, 
s, c, b, t) and six leptons (e, z/e, /i, z/^, r, z/^) which all matter is composed of. The 
forces between the fundamental particles are described by the twelve gauge bosons. These 
particles mediate three forces; the strong-, weak- and the electromagnetic force. Note that 
the fourth fundamental force (out of the four known) gravity is not included. The weak 
force which comes from the local SU{2) gauge invariance of the SM Lagrangian and the 
EM force which originates from the local U{1) gauge symmetry is unified to the electroweak 
force (described by QED) which is mediated through the four gauge bosons {W^ , Z^ , '-/) . 
This unified force is based upon the U{1) x SU{2) symmetry. The strong force comes from 
the local SU{3) gauge symmetry of the Lagrangian. Due to the non-Abelian^ nature of the 
local SU{3) gauge group, the eight gauge bosons (the gluons) that mediate the strong force, 
carry "colour" charges such that the QCD Lagrangian incorporates gluon self interaction 
involving vertices with three and four gluons. This fact complicates things a great deal. 

2.1 Interpretation of the Lagrangians 

The theory of the SM is written in terms of Lagrangians, which describe the dynamics of 
a system in a fundamental and well known way. They are scalars in every relevant space 
and invariant under transformations. In the full field theory, the Lagrangians can be used 
in a variety of ways. 50 years ago Feynman created a simple method which reduces and 
simplifies amplitude calculations done in relativistic quantum field theory. This was an 
important discovery because the amplitude is of great importance since it incorporates the 
physical picture of a system. Feynman's method is based on diagrams (called Feynman 
diagrams) which can be used to calculate the amplitude in an efficient way, using so-called 
Feynman rules. The rules enforce conservation of energy and three-momentum at each 
vertex, so that the four-momentum for the whole vertex is conserved. 

Let us investigate an easy example where, 7r° -^ 77. The pion, using standard quantum 
field notation, is written as 

n\x) = I NAp){cipe~'''"' + aJe'P-")rf=^p (2) 

where the creation a^ and annihilation a operators, create and annhilate a pion with 
momentum p, respectively. The N^o{p) is a numerical factor determined from the Feynman 

^Whenever the order of transformations matters, they are caUed non-Abehan transformations. The 
non-zero value of a commutator including the Pauli matrices that generate the local SU(2) gauge theory 
is an example of this. 



rules. The photon field can also be described in a simular way, but since it is a vector and 
not a scalar it needs an extra index. A polarization vector e with four components is 
introduced. The photon field is then represented by 



A^{x) 



NA{q){e^,ake 



-ik-x 



* ^t ^ik-x 



+ e^«fce 



)d\ 



(3) 



where k is proportional to the momentum q of the created photons. For the -k^ -^ 77 
where one pion is annihilated and two photons are created, the corresponding terms are 
picked out and used in the Lagrangian to get the Feynman amplitude for the vertex. 



2.2 Renormalization 

Start by considering lowest order ("tree level") Feynman diagrams, where A + A 
Assume that the time axis goes from the bottom to the top. 

B B 

C 



B + B. 





The diagram has two vertices, so the matrix element M. is proportional to the coupling 
constant squared {g"^). This diagram is not the only one that describes this process, e.g 
the four-vertex diagrams below also contribute to the process. 

B B^ \B B. 

C ^ \. A 

B 






C 



A^ -^ A^ 



Let us investigate the bottom right diagram containing a loop integral where the four 
dimensional element could be written as df'q = q^dqdQ, where dQ represents the angular 
part. At large q the integrand is essentially 1/g^, so the loop integral has the form of: 



rq'^dq = Ing 



CXD 



(4) 



This integral is logarithmic divergent for large q, which was a problem for a long time. 
Now there exist systematic methods to get rid of the "problem". 

The procedure goes as follows: 
(1) Regularize the integral: Done by modifying the theory, in such a way that it remains 
finite and well-defined. This can be done in a couple of different ways which all lead to the 



same physical result. In this thesis we will use two different methods, the cut-off procedure 
and dimensional regularization. The first method mentioned, is the oldest one of the two 
and it is often used on simpler cases. The other procedure, dimensional regularization, is 
a powerful method which is based on modifying the dimensionality of the integrals so that 
they become finite. 

(2) Renormalize: Where the physical parameters of interest change from the ones intro- 
duced by original Feynman rules, to "renormalized" and finite parameters containing extra 
factors. 

(3) Eliminate the regularizing parameters: One revert from the regularized theory back 
to QED, i.e. the original infinities of QED relate the original (bare) particles with the 
physical particles. This means that the physical observations of the theory expressed in 
terms of mass and charge are finite as one restores QED. 

We will now demonstrate the cut-off procedure to give a feeling about the whole topic. 
So, the first step is to regularize the integral, where a suitable cut-off procedure is used that 
removes the infinities for the moment. This renders the integral to a finite value without 
spoiling the Lorentz invariance. The cut-off is introduced as a "fudge factor" that goes 
to 1 as the cut-off —>■ oo. The nice thing is that the integral can then be calculated and 
separated into two parts: a finite part, independent of the cut-off, and a term related to the 
logarithm of the cut-off, which diverges when the cut-off goes to infinity. Now something 
extraordinary happens, all the divergent cut-off dependent terms appear together with the 
masses and coupling constants. This means that the physical masses and coupling constants 
are not the m's and g's from the original Feynman rules, instead they are considered to be 
renormalized ones, i.e 

rriphysicai = m + 5m Qphysicai = g + Sg (5) 

In the infinite cut-off limit the delta factors are infinite, but this does not matter, since they 
are not the measured ones. The measured niphysicai and gphysicai are finite, so the m and g 
in Eq. © also diverge, and makes the right side contribution finite. The remaining finite 
(cut-off independent) terms from the loop integrals, also suffer from the modifications of m 
and g. But the change in these finite effective masses and coupling constants also depend 
on the momentum of the particles involved, and therefore they are called running masses 
and running coupling constants. If all the infinities from higher order diagrams behave in 
the way described above, we say that the theory is renormalizable. 30 years ago 't Hooft 
showed that all gauge theories, including QCD and QED, are renormalizable jH]- 

Another interesting way of viewing renormalization is with the help of a medium. 
From quantum physics we know that every particle moving through a medium consisting 
of quantum fluctuations which effect values i.e. charge and mass. To be able to ignore 
the medium the values of the particles parameters have to be changed to scale-dependent 
values, and you say that the parameters have been renormalized by the medium. 



2.3 Symmetries 

Symmetries are very important in particle physics. Whenever there is an invariance there 
is a corresponding conserved quantity. That is exactly what Noether's theorem says: 
Symmetry <^ conserved quantity. These conserved quantities and the associated currents 
create, with the help of group theory the building blocks of the standard model. The SM is 
full of different symmetries. It is built up out of the three local gauge symmetries SU{3)c x 
SU{2)l X f/(l)y, some continuous global symmetries and some discrete symmetries. 

The full QCD Lagrangian with the SU{3)c gauge symmetry can be approximated by a 
light-flavor Lagrangian, because the scale of interest is ~ IGeV. So the c, b, and t quark 
will not be considered, since they are created at a scale over ~ IGeV. The light-flavor 
QCD Lagrangian can be written as 

^QCD = 5Z ^^lIPql + iqRipqR - mg{qLqR + ([nqL) (6) 

q=u,d,s 

where qi^R = |(1 ± 75)^', Ip = i'j'^id^ — iG^) and G^ the gluon field. The q corresponds to 
a particular quark state, and its wave function is a product of factors, 

space \ / spin \ / U(l) \ / SU(2) \ / colour \ 
factor J y factor J y factor J y factor J y factor J 

where each factor represents some labels, coordinates, or indices. If the masses in the 
approximated light-flavor Lagrangian are put equal, the Lagrangian gets an extra global 
symmetry. The symmetry is an S'f/(3) flavor symmetry (rotation in quark flavor space). 
In the chiral limit the quark masses mg=0, which implies that S'f/(3) -^ SU{3)l x SU{3)fi 
since the left and right handed quarks decouple from each other. So Cqq% is invariant 
under SU{3)l x SU{3)r x U{1)a x U{1)v = J- If we define 

qv 

qs 

then the transformations associated with J symmetry correspond to 



SU{3)v- 


--L+R 


^ _^ g-je.A/2^ 


U{1)a 




m -^ e-''^''-^ 


Uil)v 




m -^ e-''^) 


SU{3)l 




^L -^ Ql'^l 


SUi3)n 




^K -^ gn'^R 


SUi3)L 


X SU{3)r 


^L + ^i? ^ Ql'^l + Qr'^r 



gL,ReSUi3)L,n (9) 

where A is a vector with the Gell-Mann matrices as its elements. These transformations 
are later on used to derive the structure of the effective Lagrangians. 
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A symmetry is usually implemented in one of the two following ways: The Weyl mode 
and the Goldstone mode. In the Weyl mode the Lagrangian and the vacuum are both 
invariant under a set of symmetry transformations generated by Q, i.e for the vacuum 

e''^ I 0) =1 0) ^ Q I 0) = (10) 

where Q is a generator of the symmetry group under consideration. The vacuum | 0) is 
unique and defined as the state of the system which 

(0 I /J I 0) = min. (11) 

So the generators Q generate a symmetry i.e. when [Q, H\ = 0, and | a) and | a') belong 
to the same multiplet then H \ a) = Ea \ a) implies that H \ a') = Ea \ a'). This means 
that a and a' are degenerate states. 

In the Goldstone mode the Lagrangian is invariant under the symmetry transformations 
generated by Q, but Q \ 0) j^ 0. The consequence of this will be discussed in detail in the 
next chapters. 

2.3.1 Spontaneous symmetry breaking 

In the following we consider an easy example that describes the theory behind spontaneous 
symmetry breaking. This happens when the symmetry of the Lagrangian is not shared by 
the groundstate. We will develop this idea and consider the linear sigma- and the non-linear 
sigma-model, where an effective Lagrangian is derived from the symmetry properties of the 
theory. Later on we use the same symmetry arguments for the QCD and QFD Lagrangians 
to determine their form. 

Start by consider a Lagrangian with a discrete reflection symmetry, -^ —0, i.e the 0^ 
theory Lagrangian 

^ = ^(5.0r-^/^V-^0^ (12) 

where is a spin-zero scalar fleld, A and n are constants. The fleld 0o is the classical 
minimum of the potential 

V{<P) = 1/iV^ + ^<P' (13) 

which has two minima for /i^ < 0, 



The constant v is the vacuum expectation value of (p. Now expand around one of the two 
minima, 

0(x) = V + a{x) (15) 



V 



(^f-SWf-VX 







An "i" dim. hyperplane 



Figure 1: Shows the "Mexican hat potential" in Eq. (fTHj) where the potential energy l^ is a 
function of 0* and the plane parallel to V is thought of as an "i" dimensional hyperplane. 



and rewrite C in terms oiv and a. Dropping the constant terms, the following C is obtained 

(16) 



^ = ^(5.-f-^(-2^V 



A/i2 3 .V _4 



A 

4! -^"-4!^" 



which describes a simple scalar field of mass a/— 2/1^, with a^ and a^ interactions. The 
reflection symmetry in Eq. p2|) is no longer present and you say that it is hidden inside 
the Lagrangian. 

2.3.2 The linear sigma model 

A generalization of the discrete symmetry breaking above is the linear sigma model, which 
has a broken continuous symmetry. This is a phenomenological model for light mesons at 
low energy. The model consist out of N fields and has a 0^ interaction that is invariant 
under rotation of the N fields. The Lagrangian is 



C 



1 



1 



{d,^r - ;/^'(0T - ^m'Y 



X 

I' 



(17) 



where i=l, ..., N and is summed over. The last two terms form the potential and are often 
written as 



lA^r + ^mr- 



v{(i)') = ^u'i(j)'r + -ii(j)TY. (1^ 

The Lagrangian in Eq. ()17|) is invariant under the symmetry transformation 

(f)' -^ R'^(iy> (19) 

for an A^ X A^ orthogonal matrix R which imply that C has an 0{N) symmetry. The 
potential V has a local maximum at 0* = for /i^ < 0. The minimum (ground state) can 
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be found through ^ and thought of as a multidimensional circle. Excitations in the N:th 
direction (radial direction) require a shift away from the minimum, whilst an excitation 
along the circle, where the potential is flat, corresponds to the other A^ — 1 directions. Here 
an intuitive physical picture can be seen, and as will be seen later, the excitation in the 
radial direction corresponds to a massive particle and an excitation along the fiat potential 
gives rise to massless particles. The chosen minimum of the potential, (po, can be written 
in the following way: 



X 



(Plix) 





V 



1,. 

N 



AT-i 



(20) 



When a perturbation expansion is performed, the fluctuations around 0o have to be con- 
sidered instead of the fluctuations around 6 = 0. This means that 



(f)\x) = {ip^{x),v + a{x)) 



k 



l,...,iV 



1 



(21) 



where ip^{x) are the light meson fields and v the vacuum expectation value. The original 
0{N) symmetry is then hidden (broken) and the remaining symmetry is 0{N — 1), which 
rotate the ip''{x) among themselves. Now Eq. ()17j) can be rewritten by replacing the old 
0* with the new one in Eq. ()2H1 . The following result is then obtained 



c 



\{d,0'f + \{d,af -\{-2^^^)a^ 



-\lj?a 



2^3 



-V 

4 



A 



(^ 



fcN2~2 



((V^ 



k\2\2 



f22) 



where A^ — 1 massless meson fields and one massive sigma field have appeared. 
If N=4 ^ an 0(4) symmetry is achieved, which means that 



/ ^ \ 



l-'\ 



TTi 



TTn 



o 



\-u 



( o \ 

vr2 



(23) 



where O is an orthogonal 4x4 matrix that rotate the matrix 0. This means that the four 
primed fields are linear combinations of the unprimed fields, and that the fields transform 
linearly. If o" = f + a is considered, where the sigma field is a scalar field and has a 
non-vanishing vacuum expectation value, then Eq. (j22p can be rewritten as 



/i^ 



A 



?2\2 



Using a representation called the exponential representation 

S = a/o^ + (^ — t> , U = exp{iT ■ (p' /v) . 



(24) 



(25) 



'^N=4 is useful because 0(4) = SU{2) x SU{2) which means that the exponential representation with 
the Pauli matrices can be used. 



The resulting Lagrangian can be written in the form 

C = -{{d^Sf - 2fi^S^) - \vS^ - -S^ + i!i±^(5^[/9'^f/t) (26) 

where (...) is the trace in flavor space. Since we have done the switch — >■ t/, using 
local group arguments^, the linear transformation is not valid anymore. The matrix U is 
invariant under SU{2)l x SU{2)ji, but this is a non-linear transformation, which can be 
seen explicitly when we Taylor expand 

f/^l + lr./-i^ + ... (27) 

and use that U -^ giUg]^ =^ f ■ 0' ^ qlt ■ 0'g]i + ■ ■ ■ . So the Lagrangian in Eq. ^I^ can 
be thought of as the non-linear representation of the linear sigma model. 

2.3.3 The non-linear sigma-model 

When physics is studied at a specific energy scale only particles that can be created inside 
this energy scale have to be considered. The heavy fields do then not turn up explicitly, 
but can be seen through virtual effects. When an effective C (see chapter 2.5) is used at 
low energies the heavy fields do not have to be included, but the virtual effects are included 
via LECs connected to only light fields. There is a theorem called the decoupling theorem 
that says that all effects from heavy fields will reveal themselves through renormalization 
of coupling constants or they can be ignored because of the heavy mass suppression. At 
low energy, Eq. ()26p can be expanded to a very interesting form. By only considering 
lowest order (/^'(/^'-interactions the only interesting term is the one without the field S, so 
all diagrams with S interactions and S propagators are neglected. The non-linear sigma 
model is achieved when S — > 0, which implies that 

£ = ^(9^f/9^[/t). (28) 

This is of great importance as we will see later on. Without external fields and EM 
interactions this term is the only term of 0{p'^). 

2.3.4 Goldstone bosons 

That massless particles show up when a continuous symmetry is broken, is a consequence 
from Goldstone's theoremPJ. The theorem says that that for every spontaneous broken 
symmetry, the theory must contain a massless particle. So the result from an 0{N)- 
symmetric theory as for the linear sigma model with A^ — 1 broken symmetries, is A^ — 1 
massless particles. The massless particles are called Goldstone bosons (GB). 

"*(/?' is determined from a + if ■ = (S + v)U. 
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As I mentioned before the Lagrangian is invariant, but Q \ 0) ^ for a number 
of generators. This imphes that in the vacuum there are operators that form states of 
particles, denoted | T["'{k)). This fact will be used later on when we show that the non- 
linear sigma model is equivalent with Cqcd- 

In the earlier chapters the symmetry of the QCD Lagrangian was discussed and in 
the chiral limit a SU{?>)l x SU{?>)r x U{1)a x U{1)v symmetry was achieved. If only 
the two lightest quarks had been considered then the QCD Lagrangian would have a 
SU{2)l X SUi2)R X U{1)a x U{l)v The first two groups corresponds to 50(4) (true 
under Lie algebra). The interesting part is that since SO {A) -^ S0{?>) under spontaneous 
symmetry breaking (where 5*0(3) corresponds to SU(2) under simular group argument as 
above)^ SU{2)i x SU{2)ji — > SU{2)y=L^R, where SU{2)y=Lj^R is a subgroup of the chiral 
symmetry. Using the same arguments as above G = SU{3)l x SU{3)fj — > SU{3)v=l+r = 
H. Since the global axial part of Cqcd in the chiral limit is broken and the fact that we are 
working in SU{3), eight ^ pseudoscalar Goldstone bosons (GB) appear as a consequence 
of Goldstone's theorem. 

The massless particles showing up form the light meson octet containing: vr^, n^, K^, 
K^, K^, 7]. The reason for this is that the light particles have the correct quantum numbers 
to be associated with the generators of the coset space G/H, which parameterize the GB. 
Since the GB are massless in this theory, their masses have to be introduced in a different 
way. It can be added to Cqcd, but this will break the perfect chiral symmetries. The 
adding of the mass terms to the QCD Lagrangian is analogous to adding a a$ term to 
the potential in Eq. p8|) . Using a perturbative ansatz, one obtains massive GB, where the 
squared masses are proportional to the breaking parameter a. Since the masses of the three 
light quarks are small compared to the breaking scale of chiral symmetries (~ IGeVy, the 
explicit adding of the mass can be treated as a small perturbation. The fact that the 
pseudoscalar-meson masses are much smaller than the typical hadronic scale reinforces 
this assumption. 

One more interesting feature of this is that the non- vanishing masses of the light pseu- 
doscalars in the "real" world, are related to the explicit symmetry breaking in QCD due to 
the light quark masses which are believed to acquire their mass from the Higgs mechanism. 



■^S'C/(2) is locally identical to the rotations in three dimensions, which form the group 5*0(3), but they 
differ globally i.e. for rotations that are not infinitesimally small. The two groups have the same Lie 
algebra, i.e. their generators have the same commutator relations 

^The number of GB is determined by the symmetry properties of the groups. The G group denote the 
symmetry group of the Lagrangian, with no generators. The H subgroup, that leaves the vacuum after 
spontaneous symmetry breaking invariant, has nn generators. Every generator that does not annihilate 
the vacuum creates one massless GB, i.e. the total number of GB is nc — nu- In our case uq = 16 and 
riH — 8, which means that eight Goldstone bosons appear. 

^Around this energy the spin-one states (p, K* , (f) can be created with the same quark content as the 
light meson octet. The spin-one particles can not be GB, since GB have spin zero if Lorentz invariance is 
a good symmetry. 
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2.4 Chiral perturbation theory 

Chiral perturbation theory (ChPT) is based upon the following two basic assumptions: 

1. The spontaneous breakdown of the SU{3)l x SU{3)ji symmetry to SU{3)v for the 
QCD Lagrangian in the chiral limit. 

2. Close to the chiral limit the mass terms of the light quarks can be treated as small 
perturbations. 

It is very useful for strong interactions in the low energy sector. As seen in previous 
chapters the symmetry of the QCD Lagrangian is spontaneously broken to SU{3)v x U{l)v 
which introduces eight massless Goldstone bosons in the form of the (vr, K, rj) octet. Since 
U{l)v corresponds to the baryon number it plays no further role in low energy meson 
physics. 

When the quark masses are added to Cqcd, the pseudoscalar mesons aquire mass 
through an explicit chiral symmetry breaking. The mass is added through a chiral invariant 
coupling of the mesons to the scalar and pseudoscalar external currents (see appendix C) 
s and p, respectively, 

X = 2Bo{s + ip) ^ 2BogR{s + ip)g{ (29) 

where gL,R^SU{3)L,R, Bq is related to the quark condensate (OImmIO) = —FQBo[l + 0{mq)], 
p is a 3x3 matrix containing eight pseudoscalar external fields related to the qj^q current. 
The scalar external currents includes the masses for the three light quarks, 

f rriu \ 

s = M+-= \ ma +■■• (30) 

\ ms J 

and are related to the qq current. To get rid of the part of the external current s which is of 
no importance to us, put it exactly equal to the mass matrix. By doing this x 'hOLXdRi so 
the chiral symmetry breaks, but the quark masses are small compared to the constituent 
mass, so this can be treated as a perturbation. And since the creation of the GB is based 
upon the SU{3)l x SU(3)r symmetry which is broken by the perturbative mass, the GB 
change to pseudo-GB. 

ChPT is a non-renormalizable theory, because of the ultraviolet divergences produced 
by the loops. For a non-renormalizable theory infinite numbers of counterterms have to 
be added to make it finite. The problem can be avoided by using effective quantum field 
theory (see effective Lagrangian chapter), which indicate that the energy region of interest 
is below a certain cut-off parameter A, and then the number of counterterms needed at 
each level of expansion is finite. At increasing expansion powers, the number of needed 
counterterms increases drastically. At very low energies, the contributions from higher 
orders are small, that is why only terms up to 0{p'^) are considered here, see Eq.(^. 
The finite parts for each power contains a finite coupling constant and it is some of those 
constants that are going to be determined in this thesis. The problem is that they cannot 
be calculated directly from the QCD and QFD Lagrangians. 
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2.5 Effective Lagrangians 

The reason why an effective Lagrangian is used is that at high energies QCD is pertur- 
bative, i.e it can be expanded in terms of a^, which is small for large energies. At low 
energies the theory is highly non-perturbative since Ug is large and we therefore need to 
use effective Lagrangian techniques. When effective Lagrangians are used instead of the 
original Lagrangians, they correspond to the same symmetries but generate different Feyn- 
man diagrams. 

Most theories in physics are an approximation to the real world. The perturbative 
treatment is a common and well known method. For the task of this thesis a perturbative 
theory is needed. The problem can be solved if a light pseudoscalar meson representation 
is used. The representation is not expanded in terms of a^, but instead in terms of quark 
masses and momenta. This is known as chiral expansion. 

For low enough energies only a few relevant degrees of freedom are needed to describe 
the theory. The non-relevant d.o.f can be integrated out and encoded into the LECs, just 
as for the non-linear sigma model. The effective low energy action that is described by a 
number of relevant and non-relevant fields, can be written in the form of external fields in 
the GB low energy theory. The earlier mentioned Lagrangian is then modified, and the 
external fields can be included in a chiral invariant way. 

2.5.1 Power counting 

Power counting is a very useful way to organize the powers in the chiral expansion. The 
ordering in the context of Feynman rules is that, derivatives generate four momenta and 
the quark mass terms have the same dimension as two derivatives. So when we are on- 
shell rri^^g^ = p^, but when we are off-shell e^ has the same counting as the four momenta 
squared. The external photon fields are of order p since they occur together with a mo- 
mentum in the covariant derivative. An illustrative example can be seen below: 



meson vertex ~ p^ 




meson propagator ~ l/p^ 



/ d'^p loop integral 



4 
~ p 




(p2)2(l/p2)2p4 _ pi 
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In the figure above, the loop diagram is of 0{p'^) and is generated by tree level diagrams 
from €2- The loop integrals generate divergences which can be canceled by adding tree 
diagrams, from the £4 terms. The counterterms needed to cancel the divergences created 
by the loop diagrams of (9(p"') comes from the LECs of the same order, which implies two 
things. First the good thing which is that the loop integral divergences are canceled order 
by order. The other thing is that the number of coefficients needed to describe the theory 
increases dramatically for increasing orders, which makes an expansion to higher orders 
almost impossible. 

For TTTT scattering in pure QCD, an infinite number of diagrams have to be considered 
since they are equally important. If instead an effective theory is considered, the diagrams 
can be divided into groups of different powers of p^. For low energy, only few powers of p^ 
have to be considered. The physical picture can then be described by a small number of 
diagrams, which is a good approximation of QCD for low energy. 

2.5.2 QCD at low energy 

At low energies quarks do not behave like free quarks, instead they behave like asymptotic 
light meson fields. The effective theory of QCD is formulated in this way. It can be 
expanded in a series of infinite numbers of derivative terms, i.e 



-eff 



C2 + Ci + C(i + --- (31) 



From the properties of the GB and the symmetry of QCD ^ , the easiest form of the leading 
order terms is 

^QCD ^ ^^Qf^^^Q^jj ^ 2BoM{U^ + U)) (32) 

where M. is the mass matrix in Eq. ()3U|). U contain the eight pseudoscalar meson fields 
and has the following properties: 

UU+ = 1, det?7 = l 

f/ = exp(^), <I> = VA,<^„ (33) 

where 



a=\ 



^0 



75 + ^ -" ^" 

$ = V2| TT- -^ + ^ K^ |. (34) 

The observant reader can see that in the chiral limit, the vacuum expectation value v in the 
non- linear sigma model in Eq. (|^ can be put equal to the pion decay coupling Fq. This 

®With this I mean by adding local contributions to the external fields the transformations associated 
with the J symmetry in Eq. implies certain invariances, which are used to derive the structure of the 
effective Lagrangians. All this follows from Gasser & Leutwyler [H]. 
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coupling constant is defined from tlie differentiation of tlie lowest order classical action 
S2 = J d'^xCfz with respect to external axial fields in the chiral limit, i.e. 

(0 I ^ I 7r+(p)) = (0 I (i^Ts^ I vr+(p)) = zv^Fop^ (35) 

We compare this with the pion decay constant which is derived as follows: The decay of a 
TT"*" into two leptons via a weak interaction can be described by the tree level diagram 

^U fi^ 

...W.t... 





Written in terms of matrix elements , this corresponds to 

{W+ I C I 71+ (p)) ~ {W+ \W+ \0){0\ d^i^u I 7r+(p)) (36) 

The pion decay constant is then defined to be the second term on the right side. When 
the amplitude is calculated the following is achieved 

iV2F^p>' = (0 I d^^^^u I 7r+(p)) (37) 

and the identification F^^ = Fq can be made. 

2.5.3 QFD at low energies 

The electromagnetic virtual interactions of order e^ between pseudo-GB to the leading 
order in the chiral expansion which was first derived by R. Dashen |10J and is described 
by the effective Lagrangian 

C^M ^ (,^Ci{Q^) + e^C{QUQU^) (38) 

which is created from the symmetry properties for an effective low energy QFD theory. To 
ensure the chiral SU{3)l x SU{3)ii symmetry in Eq. (jHH|) . the local spurious^ 

Ql -^ giQgl 

Qr ^ 9rQA (39) 



are introduced instead of the charge matrix Q. The first term in Eq. ()38|) is a constant and 
will not be considered. So the resulting Lagrangian is 

£f f ^ e^dQLUQRU^). (40) 

The two spurious, Ql and Qr, are different matrices. Since the photon couples identically 
to the left respectively right, it can not distinguish between the two Q matrices and for 
that reason we put them equal in our effective Lagrangians. 

^Spurions are dummy fields that are added to the Lagrangian to make it a singlet under the chiral 
group. 
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2.5.4 Leading order effective Lagrangian 

The effective Lagrangian of QCD including EM effects to leading order is in the mesonic 

sector 

£2 = -\f,.F^'' - i^(^^(9,A^f + \f^{D^U^D,U + xf/+ + x'U) + e^C{QUQU^) 

(41) 

where, the two first terms describe the kinetic energy and the gauge fixing terni^° for the 
virtual photons. When one takes the inverse of the two terms it will result in the photon 
propagator that we later on use to calculate the photon loops. F^i. = 5^Az/ — dyA^j, is the 
electromagnetic field strength tensor of the photon field A^ and ^ is a gauge fixing term. 
The covariant derivative in Eq. (j4ip conserves the local gauge symmetries since it includes 
the couplings to the photon field A^, the external vector v^ and axial- vector currents a^, 

D^JJ = d^- i{v^ + a^)U + iU{v^ - a^) = d^ - ir^U + iUl^ (42) 

where v^ = ^v"^ = e^QA^ + -u^, and a^ = ^a^ is a 3x3 matrix associated to the quark 
current g7^7^A'*g. v^ is associated with g7^A°g, but this current is not of interest to us. 
The Q matrix in Eq. ()41|) is the charge matrix for the u, d, and s quark, 

Q = l\ -1 I- (43) 

As we have seen above, the lowest order Lagrangian, £2, is given by the non-linear sigma 
model coupled to external fields. When calculating a process using only tree level diagrams 
of C2, one represents the results of current algebra. The tree level diagrams can be put 
together to one-loops diagrams where we have an emission and absorption of a photon in the 
pseudo-GB theory, generated by the two-point function (Fig. EI), which include the vertices 
'yP~^P~ and 'y'yP^P~ incorporated in Eq. (J4H) . A calculation of one-loop diagrams with 
£2 vertices, leads to divergences as I mentioned before. Due to Weinberg's counting these 
infinities are of 0(p^) and can be canceled through renormalization by next-to- leading 
order Lagrangian at O(p^). 

2.5.5 Next-to-leading order effective Lagrangian 

The Lagrangian in Eq. (PT|) is only the first term in the infinite chiral expansion series. 
Here we introduce the second term in the expansion. The first part of it is of order p'^ 
and for the three light quark approximation in the pure strong sector {Q = 0), it was 
first written down by Gasser and Leutwyler [H], while the EM part of order p^e^ was first 
derived by Urech [7]. 

^"^The main reason why the gauge fixing term is introduced is to avoid singularities when the photon 
mass equals zero. 
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Figure 2: The two first types of one-loop contributions from Eq. ()4H). The crosses are 
external pseudoscalar currents and the wiggly line is the virtual photon. The full lines are 
pseudoscalars. 

The £4 Lagrangian contains all allowed operators that transform linearly and are in- 
variant under the correct symmetries of QCD and QFD. It is given by: 

+L,{D,UD^U^){XU^ + Ux^) + L,{iD,UD'^U^)ixU^ + Ux^)) 
+Le{xU^ + Ux^)' + LjixU^ - Ux^)' + LsixU^U^ + Ux^Ux^) 

+H,{R,,R^'' + L,,Ln + H2{x^x) 

+e^F^Ki{D^'U^D^U){Q^) + e^F^K2{D^U^ D^U){QUQU^) 
+e^F^K^{QU^D^'UQD^,U^U + QUD^'U^QD^UU^) 

+e^F^K^{{D^'UD^U^ + D''U^D^U)Q^) 
+e^F^KQ{D''U^D^UQU^QU + D^'UD^U^QUQU^) 
+e'F^Kj{x^U + U^x){Q') + e'F^Ksix^U + U^x){QUQU^) 
WF^K.,{{x^U + f/tx + xU^ + Ux^)Q^) 
+e'F^K,o{ix^U + U^x)QU^QU + (xU^ + Ux^)QUQU^) 
+e'F^Ku{ix^U - U^x)QU^QU + {xU^ - Ux^)QUQU^) 
+e^F^K,2{UD^U^[\/^Q,Q] + U^D^U[\/fiQ,Q]) 
+e'F^K,s{V^QUVlQU^) + e'F^KuiV^QV'^Q + W^QW^^Q) (44) 

L(R) 

Vj(^)Q = d,Q - i[v, - (+)a^, Q] (45) 



where the covariant derivatives V^ are defined as 
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and the L and R field strengths are 



The couphng constants Lj, Hi and Ki are defined as 

L, = r,A + L[(/i), 

i/, = A,A + if[(/i), 

ir, = S,A + K[(/i), (47) 

where the Fj and Aj constants are determined by Gasser and Leutwyler in j3] and the 
Ej constants by Urech in (2j. A is the divergent part of the couphng constants and can 
be determined through dimensional regularization and MS regularization. The L^fi), 
Hl{ii) and Kl{fi) are, respectively, finite parts, which we are going to determine. These 
coefficients are renormalized and depend on the arbitrary renormahzation scale /i. 

3 How to determine the EM coefficients 

The main idea behind the calculations of the KiS is based on the investigation of a single 
two-point function using two different methods, and comparing them in the end. On the 
RHS we use pure ChPT to 0{p^), describing both QCD and QFD with internal photons. 
The amplitudes for the relevant Feynman diagrams on this side are generated by C^ff = 
C2 + C4 + ■ ■ ■ . On the other side, the LHS, we use external photons coupled to an effective 
Lagrangian with a QCD symmetry, cfrr = C2 + C^ + ■ ■ ■ . 

First a summary of the plan: 

• Investigate the Lagrangians £2 in Eq. pij) and £4 in Eq. (jUJ. 

• Expand the U matrix and plug in the matrices in Eq. ()34j) . Eq. ()42|) and Eq. ()43p . 

• Consider the relevant Feynman diagrams. 

• Integrate over the momenta. 

• Calculate the amplitudes for the relevant processes. 

• Compare the different amplitudes generated by £^ r = £3 + ^a + ■ ■ ■ on the 
left hand side (LHS) with £e// = £2 + £4 + ■ ■ ■ on the right hand side (RHS). 

The main result in the chapters before is the Lagrangian in Eq. P^ of 0{p^) (where the 
e^p^ tree-level contribution from Eq. (j44j) also is of this order). By expanding U in terms 
of the meson matrix $ and plugging in all the matrices into the Lagrangian, the terms that 
build up the Feynman diagrams are achieved. To handle all the terms by hand is a very 

18 



bothersome task, since there are so many of them. Instead we use the computer program 
Form to calculate the terms. They are carefully investigated and divided up into certain 
Feynman diagrams, e.g tttt interaction with different kinds of loop diagrams. The terms 
from the Lagrangian are translated according to the Feynman rules and the amplitude for 
each diagram is calculated. This is easier said than done, since the diagrams with a loop 
contain a loop integral over every relevant four-momenta. These integrals are calculated 
with different methods, depending if it is on the RHS or the LHS. The reason for this is 
that the two sides are related to different theories describing the same two-point function 
(appendix B). 

We will now investigate the two sides separately. 

Left hand side, (LHS) Here we use a Lagrangian with a QCD symmetry including ex- 
ternal virtual photons. There exist no divergences except the ones that can be ab- 
sorbed with renormalization into niq, ag and a em- The following picture illustrates 
the process 




QCD 



where the figure describes a two-point function, with a virtual photon connected to 
it. The "bubble" generated by the Green function can be thought of as a propagator 
describing a process from one side to another. 

The vital point is that since we are dealing with virtual photons an integration 
over the whole four-momentum space has to be performed. This is going to give 
us divergences, which are absorbed as I mentioned before. The integrals over the 
momentum are divided up into two parts with the help of a cut-off (not the same 
cut-off used in the renormalization chapter) which respectively is determined by a 
different method. The two parts are separated into a long-distance (LD) and a 
short- distance (SD) contribution, with the help an euclidian cut-off A. The cut-off 
also works as a matching variable where a plateau in A indicates matching. After 
passing from Minkowski space to Euclidian space with the help of a Wick rotation^^, 
the integrals can be written in the following form 

/•oo /"A /"oo 

/ drE= drE+ / drE. (48) 

Jo Jo J A 

where te is the momentum. The first integral, the LD part, is evaluated using ChPT 
to 0{p'^) which we will explain more in a moment. The second integral, called the 
SD part, is calculated with perturbative QCD. 



^^This means that the zeroth component of the four-momenta, p*^ = ip%, where p^ is in the Minkowski 
space and p^ in the euchdian space 
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SD We evaluate the photons which have a moinentum larger than A. This is done 
through an expansion in powers of 1/A, where contributions up to next to lead- 
ing order are considered. The essence of the physical picture can be illustrated 
through a simple one-loop diagram. 




We use large photon momenta, so that the Q^ is big. This means that a gluon 
or a quark has to be included to run back with the momentum Q, otherwise 
the diagram is heavily suppressed. The suppression comes from the fact that if 
only "one-way" diagrams are considered, where the energy flows in one direction 
then the creation of the final hadron is suppressed since the quarks carry too 
much energy. 

In the calculations of the SD part we consider four types of different contribu- 
tions, each associated with one type of diagram and expanded to order 1/A^. 





(b) 



C 

o 
o 
o 
o 
o 
c 




(c) (d) 

Figure 3: These are the four contributing diagrams. The curly line is the gluon, the wiggly 
the photon, the plain line the quark and the cross marks the external current. 

In the leading order of l/Nc^"^, Fig. Ot and Fig. Etl are just a product of pure 
QCD currents and their contributions have already been calculated in ^J. The 
other two diagrams have to do with the renormalization of the scalar, pseudo- 
scalar, vector and axial-vector currents, and the contribution depend on the 



^^The \/Nc is a procedure where you assume the number of coulors Nc to be large so that you can 
expand in 1/Nc- 
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scale where the input current quark masses are renormahzed in QED. This 
contribution has to some extent been calculated in P], but in Ref. [12] by B. 
Moussallam it was pointed out that an extra contribution was needed. This 
implies a different SD result for the K^ to K12 coefficients. We will use the result 
from pp as much as possible in this text and just add the ChPT corrections of 
0{p^) which originates from Fig. Et^ andEb- 

LD For this part we want to include all the orders in the QCD chiral Lagrangian 
expansion, so that a proper and a smooth result will appear and all the polyno- 
mial divergences will disappear. This is obviously not possible since there are an 
infinite number of orders. The expansion up to £g [14] is known, but we will 
only consider C^ and £4 . To this order the following relevant Feynman 
diagrams are generated: 





m — ^ ^ ^ ^ — h 





f % — 53 153 i^ 



The top left diagram might seem strange, but it just represents two external 
fields coupled to a photon loop. Thus it consist out of one vertex and a photon 
propagator. 

Since we are only considering the leading and the next-to-leading order, prob- 
lems will appear. For the LD part these problems can be avoided by using low 
energies, i.e. a small A which implies that terms like Lf4^, Lf^, ... <^ Li^. 
This means that to a good approximation for small momenta the higher order 
terms do not contribute. 

Right hand side, (RHS) Here we have an effective Lagrangian in pure ChPT describing 
both QCD and QFD with internal photons. This means that ordinary field theory 
can be used. Since the tree level diagrams created by £2 are 

of order p^ and e^ 
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the vertices of O(p^) can be put together to one loop diagrams of order p*^, e^p^ and 
p^e^. These diagrams generate divergences of the same order, which can be removed 
by adding tree level diagrams from 

£4 with Lj of order p^ 
£4 with Ki of order e^p^ 



The Li coefficients will absorb the p'^ divergence from the loop diagram created by £2, 
and the Ki coefficients will remove the e^p^ divergence. Since we are only interested in 
the Ki coefficients, the order that we are going to consider is e^p^ . By renormalizing 
with dimensional regularization and MS which is just a method to isolate the loop 
divergence and include some specific constant terms, all the infinities are absorbed 
into the LECs (the Ki and Lj coefficients in this case) and the Green function is 
well defined. When this renormalization is performed, a renormalization scale /i will 
appear just as in |2j, this means that in the end the LECs will depend on this scale. 

The resulting amplitudes from both sides are polynomials and logarithms of momenta, 
mass and electric charge. The identifications are done particle by particle combined with 
different external currents. Since the loop integrals are hard to calculate, simplifications, 
i.e putting m^ = or p^ = on both sides helps the calculations. The one-loop dia- 
grams will generate infrared (IR)^^ divergences of the form log 3-^, log -^ for the LHS and 

2 2 

log 3-T, log -^ for the RHS. As one can see this IR divergent part will cancel on both sides. 
When both the SD and LD part on the LHS have been calculated and all the LECs 
of interest have been written in terms of A the two parts are plotted against each other, 
and compared to see at which region the result is approximately the same. In this way a 
reasonable region for the SD-LD separating scale A can be determined, and it can then be 
removed. 



4 Result 

In this section the result for all the calculations are presented: 

From the two-point function, using ChPT in the external field formalism, the LD part 
and the RHS can be calculated. After matching the amplitudes on the two sides the 
following relations are obtained: 

-^ = 11^ + ^ ("9' 



^■^Just an other name for a something that go to infinity when the momentum goes to zero. The ultra 
violet (UV) divergences that show up in the SD contribution, goes to infinity when the momentum goes 
to infinity. 

22 



Achieved through: (0 | T(7r+(x)7r+(y)) | 0), using external pseudoscalar currents in the 
chiral hniit to the 0{p^). 



Kl 



3 A^ 



8C 



^^(4ie-L3-^(2LS-i;) 



3 k^U 



(50) 



Achieved through: (0 | T{n^{x)7i^{y)) \ 0), using external pseudoscalar currents with 



p^ 



m = to the 0{p^). 



K^ — K^ 



1287r2 



A2 

(e-i)iog(-) 

/i2 



3A2 



(51) 



Achieved through: (0 | T{7T^{x)7T^{y)) \ 0), using external pseudoscalar currents with 



p" 



m.o 



to the 0{p^ 



-"■9 



1287r2 

1 
1287r2 



(^ - 1) log(^) + 
(e - 1) log(^) + 



3A2 

-^(64L^6-18L^-L^) 

3A2 



Achieved through: Eq. (jKUjl and Eq. (jKT|) . 



4ir[i + {Kl - 2KI 



327r2 



,A2 



^(2^a-^ 



(52) 



(e + 2)log(^) + ^-12^L^' 



(53) 



Achieved through: (0 | T(7r+(x)7r+(y)) | 0), using external pseudoscalar currents in the 
chiral limit to the 0{p'^). 



2KI, - {Kl - 2K 



A^ 



(54) 



32.4-^°^^f)+^^°s(f) + ^^ 

Achieved through: (0 | T{7i^{x)7i^{y)) \ 0), using external axial-vector currents in the 
chiral limit to the O(p^). 



2K[^ + K 



1 



12 



647r2 



A2 1 A2 
:i + 201og(-) + --6-^LS 



/^ 







Achieved through: Eq. (j3^ and Eq. (jMj) 



(55) 



- Kl; + 2Kl^ + AKl^ - 2{Kl - 2K 



1 



-21og( 

A^ 



,A2 



/i- 



a- 



4 



1 A2 
^log(^)) + 



327r2 

+ ^^9 + ^(^10 - 2//[)(12 - 30 



(56) 
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\{Kl + Kl) 


= -{Kl- 


-2k; 


k; + Ki 


= 




k; + Ki, 


= 





Achieved through: (0 | T(7r+(x)7r+(?/)) | 0), using external axial-vector currents in the 
chiral limit to the (9(p°). 



-2iKl + K;) (57) 

(58) 
(59) 

In the calculation of the coefficients from the charged currents in Eq. ()50p to Eq. ()56|) we 
have used the result from the neutral currents in Eq. (j37|) to Eq. ()59p . The observant reader 
has probably already noticed that some of the coefficients above are linear combinations of 
others. All together there exist 10 real combinations, the ones that are linear combinations 
are Eq. (j52j) . Eq. (j55|) and Eq. (j57p . The reason why we display the linear combinations 
is that only certain combinations of the EM LECs have been calculated before, and to 
achieve those we need to rearrange our coefficients. 

The processes that we used to calculate the coefficients above are not the only ones 
investigated, e.g. for the charge current we also investigated K~^, but the combinations 
that came out of those calculations gave the same result as for the vr"*". For the neutral 
current we investigated contributions from rj, ir^, K^ and K*^ , which all contributed to the 
combinations in Eq. ()57|) to Eq. ()59|1. 

For the SD part which we calculated with perturbative QCD we got after comparing 
to the RHS: 



C = -B-^ (60) 

Kl = ^^-:^^ (61) 



3as 


F^Bl 


Stt 


A2 


3as 


F2 

-^0 


327r 


A2 







-3«. F^ 



Stt A2 
Kl = (62) 

^ 327r A2 ^ ' 

Kl = (64) 

«£_-?£ _ Pis L^Bi 
'67rA2 Svr A2 
's LIB^ 

27r A2 

Stt a? 

r r TD' 

TT A2 



_ a^F^ as LlBl 

^' ~ 6^A^ 3^^^ ^^^^ 

^6 - 2vr A2 ^^^' 

Kl = ^^ (6^ 
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KI, 



K 



10 



-"■11 



-"^12 



-a. 


, (2Lg + Hl)Bl 


67r 
3a. 


{2LI 


A2 

+ ^20i?o' 


47r 
3a. 


{2LI 


A2 



47r 



A2 



1287r2 ^A2 
1287r2 ^A2 

-(l^iogii 

1287r2 ^ A2 



647r2 ^°^A2 



(69) 
(70) 
(71) 
(72) 



where the log contribution conies from Fig. ^. and Fig. Eb- 

To be able to plot the SD part against the LD part, all the relevant coupling constants 
have to be known and the renormalization scales have to be chosen. All this is displayed 
in Tab. 0] As one might have noticed there are large errors in most of the parameters, and 



Coefficient 


fit 10 


Oip') 


article 






lO^L^ 


0.43± 0.12 


0.38 


i 




lO^Ll 


0.73± 0.12 


1.59 


IH] 




lO'-'Ll 


-2.35± 0.12 


-2.91 


IH] 




W'Ll 








IH] 




lO'Ll 


0.97±0.11 


1.46 





Bo[GeV] 


1.63 


lO'Ll 











as 


0.3 


lO'^Uj 


-0.31± 0.14 


-0.49 





fi[GeV (ChPT scale) 


0.77 


lO'Ll 


0.60± 0.18 


1.00 





fio GeV (QCD scale) 


1.0 


lO'Ll 


5.93± 0.43 


6.0 









lO^L^o 




-5.5 







103(2L^ + H^) 




2.9± 1.0 







10^ H{ 




-4.1 






Fo[MeV] 


87.7 


81.1 


[3 





Table 1: Displays the coupling constants that are needed to calculate the KiS. The finite 
renormalized parts are at the scale Mp = 0.77 GeV. The values of the parameters to the 
left are all determined from articles. As one can see there are missing three fit 10 values. 
For these we use 0{p^) in both cases. The parameters to the right have values that we have 
chosen ourselves. Bq are on the QCD scale ~ 1 GeV and a. is assumed to be constant. 
The ChPT parameter fi that appeared when we performed the dimensional regularization 
is chosen to be /i = Mp = 0.77 GeV. 

that will of course affect our result in the end. This is the reason why we have chosen two 
different fitting methods, fit 10^^ and 0{p'^Y^, respectively, i.e. the change in the result in 
the end can be seen. This will of course affect some KJ constants more than others. In 



^''include terms up to 0{p^), and is based on that L4 and Lg is equivalent to zero. 
""^^Include terms up to Oip"^), and is based on that L4 and ig is equivalent to zero 



25 



our calculations we also fix the gauge of the photon and use Feynman gauge, which means 
that the gauge fixing term ^ is put equal to zero. 

All that remains now is to compare the LD with the SD, component by component, 
and that is done in Figs. ID to Fig. ^2 On the y-axis we have the value of the LEG and on 
the X-axis we have the cut-off A. The scale on the x-axis range from 0.3 GeV to 1.1 GeV, 
and the reason for this is that a A with a too high energy (~ 1 GeV) would disturb the 
perturbative treatment (can be seen from the figures where the higher order deviates faster 
for bigger A) and since our ChPT scale is around 1 GeV, it is not proper to exceed that 
energy. The reason why at low energies it does not work is that the perturbative expansion 
is not valid below ~ 0.3 GeV. 

The matching is done by finding the best region for LD4 and SD. This is done by 
looking at the region where the Sum4 function has the derivative closest to zero, but one 
also has to consider how the LD and SD parts behave. From Fig. |3] and Fig. El describing 
the coefficient C one can see how the next order (LD4) affects the result. If only the lowest 
order LD term were considered, the matching region would be easy to find. By adding the 
next order one can see that the matching region is not so good anymore, which means that 
the result would have a larger error. The difference between LD2 and LD4 is important, 
since in this text we are also interested in comparing the behavior of the two functions. As 
one can see there is no big difference between the fit 10 input and the 0{p'^) input, that is 
why we only include the fit 10 graphs for the remaining quantities. 

When K^ is calculated as in Fig. IHlthe L^ and the Lg were put equal to 10~^ to avoid a 
non-zero result. In the figure we see that the sum of the LD and SD behaves nicely, which 
is good. 

For Kg in Fig. [7| both the sums behave nicely. The best matching region is a little too 
high, but it is still good enough. When we compare the two sums, we see that Sum4 is a 
little better. 

For the 4i^[^ + {K^ — 2KI) combination in Fig. |H1 the Sum4 function behaves in an 
appealing way. The behavior of Sum2 is also rather good, the only thing to complain 
about is that the zero derivative area is a bit too low for Sum2 and a bit too high for 
Sum4. There is also a large difference between them. 

Next we consider Fig. il Fig. HU] and Fig. HU describing 2X1^ - {Kl - 2KI) , -K^^ + 
2Kl^ + 4:Ki2 — 2{Kl — 2KI) and 2K{^ + ^[35 respectively. These graphs do not have a 
Sum4 that behaves as we want in the region of interest, but we still pick a A in our region, 
which means that the error in these parameters are big. For the 2i^[2 — (-ft'4 — 27^3) and 
—Kl^ + 2Kl^ + 4:Ki2 — 2{Kl — 2K^) combinations one can also see that the Sum2 behaves 
much better. 
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C using fit 10 input 
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Figure 4: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for C The number after the functions represents which 0{p"') that 
are considered. The fit 10 input parameters are used. 



C using 0(p ) input 



O 




0.4 0.5 0.6 0.7 O.i 

A [GeV] 



1.1 



Figure 5: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for C The number represents which 0{p"') that are considered. The 
0{p^) input parameters are used. 
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Kg using fit 10 input 
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Figure 6: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for Kg. The number represents which (9(p") that are considered. The 
fit 10 input parameters are used. 



Kg using fit 10 input 
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Figure 7: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for Kg. The number represents which 0{p'^) that are considered. The 
fit 10 input parameters are used. 



28 



4K„+(K4-2K3) using fit 10 input 
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Figure 8: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for 4:Kli + {Kl — 2K^). The number after the functions represents 
which 0{p'^) that are considered. The fit 10 input parameters are used. 
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2Ki2-(K4-2K3) using fit 10 input 
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Figure 9: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for 2X^2 — [Kl — 2K^). The number after the functions represents 
which OIp"') that are considered. The fit 10 input parameters are used. 
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-Ki3+2Ki4+4ki2-2(K4-2K3) using fit 10 input 
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Figure 10: The long distance (LD), short distance (SD) and the sum as a function of 
the matching variable A for Kl^^ + 2Kn — 2{Kl — 2KI). The number after the functions 
represents which (9(p") that are considered. The fit 10 input parameters are used. 



2K11+K12 using fit 10 input 
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Figure 11: The long distance (LD), short distance (SD) and the sum as a function of the 
matching variable A for 2i^[^ + K12. The number after the functions represents which 
0{p"') that are considered. The fit 10 input parameters are used. 
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The result from all the figures are summerized in the Tab. |21 below: 



Table 2: The results for fit 10, "This work 1" is the value calculated in this text using fit 
10, for "This work 2" we are using O(p^). The A is value of the cut-off where the SD part 



best matches the LD part. 










Coefficient 


This work 1 


A, fit 10 

[GeV] 


This work 2 


A, 0{p') 
[GeV] 


C 


3.0 X 10-5 

[GeV] 


0.62 


2.25 X 10-5 

[GeV^] 


0.58 


K'^ (used Ll = Ll = 10-3) 


5.8 X 10-3 


0.72 


6 X 10-3 


0.71 


K^9 


-8 X 10-^ 


0.8 


-9 X 10"^ 


0.73 


AK[, + {Kl - 2KI) 


-2.3 X 10-2 


0.83 


-2.3 X 10-2 


0.71 


2K[^ - {Kl - 2KI) 


4 X 10-3 


0.65 


5 X 10-3 


0.65 


-K{^ + 2K{^ + 4i^i2 - 2{Kl - 2KI) 


1.1 X 10-2 


0.6 


1.4 X 10-2 


0.6 


2KI, + Kl^ 


-8 X 10-3 


0.8 


-8 X 10-3 


0.8 



Table 3: Displays the relevant calculated coefficients for each article. For [I] -Fq = 89Mel^ 
was used, for ^2] the improved Z value for K^ and Kiq was used. In all the calculation of 

the values we have used Ho = I GeV, My = i^ = Mp = 0.77 GeV, F^ = Fq = 92 A MeV, 
Ma = Ma, = 1.23 GeV and leading Nc order. 



Coefficient 


[\ 


[12J 


[131 


C [GeV^ 


4.2 X 10-5 


2.36 X 10-5 




Kl 


0.8 X 10-3 







Kl 


-1.3 X 10-3 


-3.62 X 10-3 




AKl, + {Kl - 27^3^ 


-5.0 X 10-3 






2KI, - {Kl - 2KI) 








-Kl^ + 2K[4 + 4i^i2 - 2{Kl - 2KI) 








2KI, + Kl^ 




-5 X 10-3 


-1.94 X 10-3 



The relevance of some of the results in Tab. El is now discussed. 

When we compare the values for C in Tab. |21with the values in the articles in Tab. El 
one can see that they are rather simular, which is good. 

The value for K^ will not be discussed because of the arbitrary chosen L[ parameters. 

Instead we turn to K^ which have a good matching region, and as we see these values 
are rather close to the one derived in pp. 

The values for 'iK^^ + {Kl — 2KI) are only calculated in PP, and they differ a lot 
compared to our values. Why is hard to say, but there are a lot of approximations that 
have been done in both cases. 
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For 2K[2 - {Kl - 2KI) and -Kl^ + 2Kl^ + AK12 - 2{Kl - 2K^) no published values 
are accessible. So their values will not be commented further. 

The last combination obtained is 2Kli + KI2, which has values that are close to the 
value calculated in |12j|. This means that our values are reasonable^^, despite the large 
uncertainty. 

In Tab. 121 one can see that the fit 10 and 0(p"^) input do not differ much, which means 
that the largest error in the K^ coefficients are not due to the uncertainties in the L[s, but 
instead due to the matching between the LD and SD parts. 

5 Conclusion 

In this thesis we have determined the renormalized and finite coefficients of the generating 
functional of ChPT up to 0{p'^) including virtual photons to one-loop. The finite part in- 
corporates 29 low energy coupling constants, and of those, 10 combinations of some the 14 
K^s have in the MS scheme been written in terms of the experimentally known constants 
LI, the cut-off A and the dimensional regularization parameter /i. The A dependences 
have been removed through the matching of the LD part to the SD part, which had a good 
match for some of the coefficients and a not so good for others. But over all the 0{p'^) 
contribution made the result a little bit better. But the improvement is still not that big as 
we hoped it would be, which implies that the result of the naive leading order estimate in 
fn\ is often good enough for a first estimate in the determination of the Ki constants. One 
could argue that a contribution of 0{p^) would improve this result, but it seems unlikely 
since higher orders tend to deviate more for As in the upper region. 

The determination of the effective low energy chiral Lagrangian is of great importance 
since it is used to determine the free quark mass difference mu — rrid, which is not very 
well known. The procedure used, is based on that the physical mK+ — rrixo mass difference 
has a pure QCD contribution due to the difference in quark masses and an electromag- 
netic contribution. Since the mass of the mesons can be determined experimentally, and 
quite accurate, then by calculating the EM contribution, the quark mass difference can 
be determined. This contribution can be calculated through different ChPT techniques, 
some better than others, often depending on region of interest. The procedure used by 
Bijnens-Prades in P] includes all orders of p^ in ChPT for the LD part, and this is done by 
using the ENJL model. Ananthanarayan-Moussallam in [^ and JH] also used a method 
where all orders of p"^ in ChPT were considered. The difference is that Moussallam used 
resonance saturation instead, but we will not explore the two methods more since it is 
beyond of the scope of this thesis. 

One of the problems calculating the next-to-leading order electromagnetic coefficients 
is that the QCD coefficients are needed, the L[s in this case. Most of these terms are not 
so well determined. In the near future, hopefully, these constants can be determined more 



""^^To say that our value is reasonable, does not mean that it is the right value, since we do not really 
know which value that is the right one. 
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accurately. 

Using the ChPT expansion on two-point functions as we did for the LD part, was not 
enough to determine all the LECs, as we were hoping. We even did the calculations for 
the four-point function for the K^ {K ), but no new relation appeared. In the future one 
should investigate the four-point functions further and see if new i^[ relations drop out. 
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Figure 12: One of the basic Feynman diagrams for a two-point Green function. 

A An easy example 



Start by considering the leading-order effective QCD Lagrangian in Eq. (jHzj) . The U field 
in Eq. ()33|) and Eq. (jM|) is expanded in the following way: 



2 



f/ = exp(-)^l+.--^ + 



f^^ = exp(-)^l-.— --^ + --- (73) 



— ) ~ 1 — i 

where $''' = $ has been used. Taking the partial derivative of U in Eq. ()73|1 yields 



d,U^ = -yd^'^ - ^{d.m + $9^$). (74) 



Putting Eq. (f7^ into the covariant derivative in Eq. (jl^ gives 



M 



i^iC/ = -f a^$ - ^(9,$$ + $9,$) - ^e[QA„ i^$ - |!] + • • • (75) 



Now an expression for the first term in the Lagrangian in Eq. (jH^ has been derived and 
the Feynman diagram in Fig. El is ready to be evaluated. 

1. Here a three vertex is investigated, where the only relevant terms from Eq. ()32p are 
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the ones with one photon field A^, one $ term and one 9^$ term. So 

— 2p?/1'^ 
tr{D,UD^^U) = ■■■ = -—^tr{d,^Q, $]) + ... = ..• = 

-^(-9^7r+7r- - d^K+K~ + 8^11-71+ + d^K' K+) + ■■■ . (76) 



i^d^ 



This is interesting, since only tt and K mesons showed up. Using ordinary Feynman 
rules and multiplying in the constants 

A>^d,T,-^+ - A^S^TT + TT- ^ l{pl + ql) + tpl 



F^ Aei 



-± — {ql + 2pl) = ei{ql + 2pl) (77) 



^e^(gf + 2p^). (78) 

2. The same as above, except that fi ^ u. 

3. Here one external field couples to an interior. The second term in Eq. (j32j) contain 
the field couplings needed. The 7r°, K^ and rjs is not considered since they did not 
show up in the three vertex. The resulting terms are then 

trixU^ + Ux^) = tV2FoBo{P^+n- + Pk+t^') + ■■■ (79) 

where 

iV2FoBoP^+7i- -^ iV2FoBo 
iV2FoBoPK+K- -^ iV2FoBo. (80) 

4. The same result as above 

5. In the propagator the derivative term squared and the mass term are of interest. So 



tr{^{d,<l>y)+tr{xU^ + Ux^: 



'Fo 



45o 



-^{2d^TT+d^7T- + 2d^K+d^'K-) ^{2rmi+Ti- + mK+R-) + ■■■ (81) 

where m = (?ti„ + md)/2. The complex fields above propagate in the following way 

p; — 2Bom 

d^K+d^K- - Boifh + ms)K^K~ -^ -|— ^. (82) 

p^-5o(m + m,) 
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6. Achieve exactly the same terms as above 

7. Here the resulting terms are 

i 



{Pn + qny - 2Bom 
i 

{Pk + QkY - Bo{m + rus) ' 



(83) 



8. This is the photon propagator, and only terms including —^F^yF^^" — ^.^_^A d^A^y 
are considered. The terms propagate according to 

Adding all the parts for the Feynman diagram gives the following amplitudes for the tt 
and K case, respectively: 



■ [(?-(g^+2p^)] 



(85) 



(27r)4 (p2 _ 25om)2((p^ + q^y - 2Bom)q^ 

dS -2FiBle\{qK + 2pKf - e '^^'^^^+^^^^l^ ) 

(27r)4 {pI - Bo{m + m,)y{{pK + qx? - B^{m + m,))g2 



In the next step we calculate the one-loop integrals, where we integrate from zero to infinity. 
After transforming from Minkowski space to Euclidian space [p^ = {p^Y ~ {pf' —^ ~p\ ~ 
{p^Y + (p)^) "^6 can introduce a cut-off A and split the integral into two parts. The first 
part from zero to A gives us the following integral where the q independent terms are 
included in B 

io (27r)4 {{p^ + q^y-2B,m)q^ ' ^ ^ 

For the photons with a momentum greater than A you have to use other methods to 
calculate their contribution (see section 3). The important part is that in the end we want 
to compare the p^e^ contribution from the effective QCD Lagrangian in Eq. (jH^ with the 
p^e^ part of Eq. (jH)). This is done by putting e.g niq = or p^ = on both sides, i.e. the 
integrals are greatly simplified. 

B Generating functionals 

The two following sections are based on Ref. [T^ . 
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In this thesis Green functions in the form of two- respective four- point functions are 
used. They both originate from the matrix elements describing the physical process and 
the general matrix element is connected to a n-point function 

(0 I T(0(x.) ■ ■ ■ 0(x,)) I 0) = z«^-i!i^i^^^ = G(")(xi, ..., x„) (87) 

where j(x) is the classical source field used to probe the theory, and it can be thought of 
as an external field. A physical picture of the above is shown in Fig. ^1 The generating 



Xi 

G'(")(xi,...,xO = ^2 



Figure 13: The figure shows n external fields coupled to some process, 
functional W is associated with the sources in the following way 

W[j]= /"[rf0]e^/'^M^(^.9'/')-^-<^). 




and one should think of C in Eq. (|Hg|) as describing physical processes inside the big bubble 
in Fig. ^1 and the external fields as creating and annihilating particles on the edge of the 
bubble. So the bubble contain all relevant Feynman diagrams, e.g the ones in Fig. El 

The two-point function is often written as a Fourier transform into momenta of the 
matrix elements with only two sources (four sources for a four point function). In the path 
integral formalism the two-point function can be written as 

n„=)^,/Ae-(0|T(P(0)pt,.))|0) (89) 

where P(x) is a pseudoscalar source in our case. 

C External fields 

Here the idea about the external field formulation using generating functionals is presented. 
The external fields are a method that we use to be able to go from the generating functional 
to the Green function, and they can be thought of as the generalized j current of appendix 
B. One can also use the analogy where a particle is accelerated into a target particle and 
we investigate what comes out. The particle in this case is our external field and the target 
is the Green function. 

The use of external fields has some advantages, i.e. it is independent of the G/H 
parameterization, allows for a well defined off-shell amplitude and it is obviously invariant 
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under chiral transformations. One more advantage is that the connection between QCD 
and the local properties become clearer. 

The external fields are added to the Lagrangian and the generating functional can be 
written as 

where the approximation is valid at low energies as we saw in the case of the non-linear 
sigma model. As seen in Eq. (jUUI) . the external currents contain quark currents and when 
the Lagrangians are differentiated with respect to the external field variables the pure 
quark current show up. The EM part could also be added to the Lagrangian above so that 
EM currents show up. 
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